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Abstract 

In a thermodynamical model of cosmological FLRW event horizons for 
the dark energy (DE), we consider a logarithmic corrective term in the 
entropy to which corresponds a new term in the DE density. This model 
of A(t) in an interacting two-component cosmology with cold dark matter 
(DM) as second component leads to a system of coupled equations, yield- 
ing, after numerical resolution, the evolutions of A(i), the Hubble H(t), 
vacuum density QA(t), deceleration q(t) and statefmder R(t) and S(t) pa- 
rameters. Its results, compatible with an initial inflation and the current 
observations of the so-called "concordance model" , predict a graceful exit 
of early inflation and the present acceleration and solve in the same time 
the age and coincidence problems. Moreover they account for the low-/ 
CMBR power spectrum suppression. 

keywords : Dark energy theory, Snla constraints, CMBR (low mul- 
tipoles) theory. 

1 Introduction 

The observations of SNe of type la [1] [2] , showing that the universe has probably 
recently entered on a phase of acceleration, raised up a renewal of interest in 
the cosmological constant (CC) problem [3] [4]. 

Since then an increasing number of observational evidences for such an ac- 
celeration has been settled. This acceleration together with other cosmological 
indicators such as the flatness (k = 0) from the CMB anisotropics observed by 
Boomerang, Proxima and WMAP [5], the present values Hq by the HST key 
project [6], £lbo and fl m o by observations of the light elements, SBBN and clus- 
ter baryon fraction from X-ray emission [7]), led to the so-called "concordance 
model" [8] based on the ACDM model. This model is the simplest improvement 
of the old standard cosmological model (including inflation) . More precisely, the 
simplest theoretical candidate explaining this accelerated expansion is a cosmo- 
logical constant A, or some component, called dark energy (DE): any term in 
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the energy-momentum tensor which can produce the same negative pressure, 
such as a scalar field [9], called quintessence [10], or a bulk viscosity [11]. 

However, the re-introduction of the CC A leads to new intriguing puzzles. 
Particularly, we have to face the old CC problem: the discrepancy between the 
predictions of the effective Quantum fields theory concerning the large vacuum 
energy density and the tiny upper limits observed in astrophysics for the CC. 
Besides, it gives rise to a new issue, the so-called "coincidence problem" , where 
one has to explain why the CC starts recently to dominate the matter. Assuming 
A = constant, it can only be solved by an unsatisfying fine-tuning. Finally, the 
introduction of A needs also to be conciliated with the inflation in the early 
universe. 

In these conditions, we have to improve the AC DM model in order to solve 
these puzzling issues. 

Therefore, we propose a A(t)CDM model, namely the simplest modification 
of the ACDM model which allows the best account of the theoretical expecta- 
tions and the observed data. Some "ansatze" were already proposed to solve 
the CC problem before the discovery of the present universe acceleration (for 
a review, see [12]). By contrast, the present model, namely the event horizon 
Thermodynamics (ehT) model, is derived from a thermodynamical approach 
and not based on an ad hoc hypothesis for the law A(t). Such a phenomeno- 
logical model is an indispensable preliminary to outline the main laws that any 
further treatment at the statistical or/and quantum level will have to reproduce, 
as for example the Stefan's law has guided the Planck law's advent for the black 
body radiation. We thus proposed a model with a variable vacuum energy den- 
sity [13] [14], based on the generalization of the Thermodynamics of the black 
hole event horizons (BH e.h.)[15] to the Fricdmann-Lcmaitre- Robertson- Walker 
event horizons (FLRW e.h.)[16]. This model allows to reconstitute the entire 
evolution from the early inflation to late acceleration (for another approach, see 
also [17]). 

In this article, we extend this phenomenological model by introducing a log- 
arithmic correction in the A-entropy. Such a corrective term in the entropy 
was recently computed for the entropy of any thermodynamical system and 
applied to Black Hole by S.Das et al. [18]. This correction due to statistical 
fluctuations around an equilibrium state has an universal character. As the cos- 
mological FLRW c.h. presents notable thermodynamical analogies with the BH 
e.h. (albeit some differences [19]), we examine the effects of the introduction of 
this new term in the A(t)CDM model. It appears that the derived expression 
for A(r), with r the radius of the event horizon (e.h.), can then be naturally 
interpreted as a renormalisation relation and allows to describe the main ob- 
served features of the evolution of the universe. In this scheme, the obtained 
universe is initially submitted to a local phase of de Sitter inflation (Oa(0) = 1). 
The vacuum density parameter begins to decrease and there is a graceful exit 
of inflation. The parameter passes through a small minimum, then grows 
again, which allows to solve the "coincidence problem" . It tends asymptotically 
to a de Sitter universe. At any instant t, H(t) and A(t) are decreasing functions. 

The paper is organized as follows. In Section 2, we set the basis of the 
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general relativistic framework in which our model for A(t) will be developed. In 
Section 3, we present the ehT model for the vacuum (or DE) energy density. 
Then in Section 4, we examine the influence of the entropy's correction term in 
the vacuum energy density and its relevance in the renormalisation approach. In 
Section 5, we deduce the equations of evolution for H(t) and r(t), and present in 
section 6, their numerical resolution, in particular, the graphic of Cl\(t). Section 
7 is devoted to a brief conclusive summary. 

2 An interacting two-component model 



Let us briefly recall the minimal assumptions underlying our model. 

In order to set our notations, let us first write the Einstein equations: 

R 

RaP - -^9afi+ A 5a/3 = X T a0- (1) 

with x — and A a positive constant (the GR bare cosmological constant). 
Only for this choice of coefficients in Eq. (1), the Bianchi identity V pT al3 = 
is satisfied [20] [21]. 

Assuming a perfect fluid of energy density pr and pressure Pt, the energy- 
momentum tensor is given by T a p = pTU a up — PrA Qj g with A a p :— g a/ 3 — u a up, 
and the index "T" means "Total". 

Following [13] and [14], the perfect fluid is supposed to admit two compo- 
nents: pt = Pi + P2 and Pt = Pi + Pi- 

The first one corresponds to vacuum, i.e., it obeys the "pseudo-barotropic" 
equation of state : P\ = upi with w = — 1. 

The second component obeys a barotropic equation of state given by P2 = 
(7 — l)p2- In section 6, we consider dust (7 = 1). 

In accordance with the WMAP observational results of the CMB's anisotropies 
[5], the metrics is taken to be the spatially flat (k = 0) FLRW metrics. With 
these assumptions, Equation (1) reduces to 

R 

R af) ~ -^9afi = kg a /3 + X T 2a0 (2) 

with T 2a/ 3 := p 2 u a p - P 2 A a/3 and 

A := A vac - A . (3) 

As A vac :— xPi with p\ variable, A varies. In the following, we refer to A as 
the " variable cosmological constant". It can be considered as an "effective" 
CC. The corresponding dark energy (DE)-component has the energy density : 

We do not make any further assumption on a separated energy conservation 
for each component. In the old standard cosmological model, this assumption is 
made only for each era of "domination" of one component (e.g. "era dominated 
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by the radiation" , "era dominated by the matter" , today "era dominated by 
the DE" , or in the very early universe "era dominated by the inflation (or the 
inflaton)" ) . If such an assumption can seem reasonable for these domination eras 
(and so appears more as an approximation), it is not the case for the transient 
periods, as for instance the very recent epoch of transition from dark matter 
(DM) to DE domination. 

Anyway, the complete energy conservation due to the Bianchi identity is 
more general than the separated conservations, which necessitate some sup- 
plementary assumption. For example, in the A(t) and quintessence models, a 
separated conservation law for each component implies a "variable" equation of 
state. By contrast, the present model assumes the equations of state in a ther- 
modynamical sense, i.e. given once for all, with an interaction between the two 
components. This kind of approach with components interaction has already 
been considered (see e.g. [22] [23]). 

Let us sketch the rough scenario of the universe evolution. Initially, setting 
Pi = PMax and p2 = P2min — 0, the universe undergoes a de Sitter inflation, 
with Am = A-vacM — A > 0. Then p\ decreases or equivalently, A decreases 
with the universe expansion. In the same time, P2 is increasing. Finally, pi 
tends to A or equivalently A tends to A m ; n ~ , which is again a de Sitter 
evolution. Today, we observe a very small value for A = xPi ~ with a 
very large vacuum energy density p\ given by quantum field theory (today pi 
is of the order of Ao). For a reasonable (Planck) largest cut-off, we obtain 
= 1 — — ~ 10~ 120 . Assuming a variable cosmological "constant" solution 

XPl XPl b b 

(as considered here), this scheme appears as an hardly intelligible fine-tuning. 
It is the so-called "old problem" of the cosmological constant (CC). Besides, the 
recent domination of the vacuum energy by dust constitutes the "new puzzle" of 
the CC (or "coincidence problem"). At the beginning, p 2 represents the energy 
density of the radiation and the ultrarelativistic matter. For sake of simplicity, 
we shall only consider in the following (section 6) the non relativistic matter 

(7 = 1)- ^ 

We propose a phenomenological model of this DE or, equivalently, this vari- 
able cosmological "constant" (or this vacuum energy density) able to produce 
the main properties of this scenario. More particularly, the two problems of the 
CC are solved in the framework of our model where A runs from a large initial 
value Am (a huge inflation in the very early universe) to a weak final constant 
A m which yet starts again to dominate, leading to "the coincidence problem". 
In the two cases, the same DE, with a density A is used to explain the universe 
acceleration. This DE is negligible or dominates, depending on the considered 
period, and its density decreases all along the universe expansion. 
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3 Thermodynamical model of the cosmological 
FLRW e.h. 

The variable DE density A = XPA is the first component of the two perfect 
fluids constitutive of the universe. In the de Sitter space-time, this cosmological 
"constant", or vacuum energy density, obeys the equations of state [15] [24]: 

A + X Pa = (4) 
Ik 2 

A = %Tl (5) 
c^h 

where Pa and T\ represent its (negative) pressure and its temperature respec- 
tively, h is the Planck constant, c the light velocity and k the Botlzmann con- 
stant. Each thermodynamical variable is linked to the proper radius r of the 
event horizon at a given instant by the relations 

<«> 

We can also introduce the horizons number density n\ defined as usual in Ther- 
modynamics by the inverse of the specific volume of the horizon [13] 

Introducing the expression (6) of the temperature in the preceding relation, we 
obtain a third equation of state: 

Equation (9) together with Eq. (4) and (5) define a local equilibrium state for 
the DE component in the de Sitter spacetime. Let us now consider the same DE 
component in the spatially flat FLRW space-time. The three thermodynamical 
equations of state (4), (5) and (9) have to remain valid in FLRW spacetime 
because the thermodynamical equations of state of any actual component (e.g. 
P = |, or the Stefan law p — o~T 4 } for the equilibrium radiation) are always 
supposed to be independent of the chosen spacetime. Besides, the definition (8) 
of «a remains also valid in the (spatially flat) FLRW spacetime. Consequently, 
from (8) and (9), we find that this DE component is associated to the FLRW 
e.h. via the relation (7). Using (7), (8) and (9), we obtain the expression (6) 
for T\ and (5) for A in the FLRW spacetime too. 

Relation (7) has also been independently derived from holographic principle 
[25] . Albeit the choice of r is not necessarily a priori given in the holographic 
approach (e.g. the Hubble horizon H" 1 can be chosen instead [26] [27], or the 
particle horizon r p as well [28]), it is the unique choice in this framework that 



5 



leads to an admissible equation of state [29] and which is strongly supported by 
the Snla observational results [30] [14]. 

Differentiating the definition of r with respect to the time (e.g. see eqs.(3.4) 
and (3.6) in [13] or (9) and (10) in [14])) yields the relation 

H- T -= C - (10) 
r r 

where H is the Hubble parameter and the dot means the time derivative. 

In order to set a complete local thermodynamics, we need to determine the 
entropy density ns A where s A is the specific entropy. To reach this goal, we use 
the Gibbs equation which at the specific level is given by: 

T\ds A = de A + P A dv A (11) 

where £a = — — is the specific energy, v A = — the specific volume. With Eq. 

n-AX ™ A 
(4), we deduce 

X n A T A ds A = dA. (12) 

4 Fluctuations and renormalization 

In a recent article, a logarithmic corrective term in the expression of the entropy 
of any thermodynamical system was derived by S.Das et al. [18] . This correction 
due to small thermal fluctuations around an equilibrium reads for BH entropy: 

s = nr 2 — a' In nr 2 (13) 

where a' is a constant. Introducing the same correction in the cosmological 
eh model, the expression of the specific entropy becomes 

s A = -4(-^-)(7rr 2 + 27raln7rr 2 ) + C, (14) 
L p 

where Lp is the Planck length (defined by L 2 P := xhc = Alp, with tp := S£ in 
the I.S. units), a is a constant of dimension [r 2 ] and C an arbitrary constant. 
As usual, the sign of the entropy is negative for the cosmological horizon (e.g. 
[16] [13]). As noticed by S.Das et al. [18] and Major and Setter [31], when there 
is not conservation of the number of particles in an open sytem, the sign "-" 
before the a-term has to be replaced by a sign "+" . This is precisely the case 
here for our A-componcnt, which is not separately conserved. 
Let us split the constant C into 

C = 47r-|-(2aln7rr^ + r 2 ), (15) 

where r 2 is associated to the corrective term and to the expression of s A 
without correction. The expression of s A becomes 

SA = 8 ^ a( ^27T _ln ^ ) ' (16) 
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The constant factor a := r\ has the dimensions of a surface. Let us remark 

the omission of the Bckcnstein factor in most papers by choice of the Planck 

( P 

units (h = c = k = G = 1). The expressions of r 2 and r 3 are constrained by 
setting 

, Vr > r(0) (17) 

where r(0) is the initial event horizon (see Section 6). The simplest choice for 
r 2 and r 3 is given by r 2 = r 3 = tm , where tm is the maximum radius (possibly 
infinite). Hence, sa is always positive with s\ = when r = tm- For r = tm — e 

with e <C rM, we deduce r§ — r 2 ~ 2er M and 2r 2 ln( M ^ rM ) ~ — 2r 2 ^-. 
A corrective term is only obtained if n <C rAf as will be verified in Section 6. 
Finally, for r = r m , we have (r 2 M - r 2 m ) - 2r\ ln((^-) 2 ) ~ r 2 M + 4r 2 ln(^). For 
fi *C I'm , the ln-term is a corrective term and tends to s — > % = ^rr-rij. It 
follows that ri has not necessarily to be very small with respect to r but only 
with respect to r M . 

Let us examine the effects of such a new term in sa (14) on the expression 
of A . 

Differentiating Eq.(14) for sa with respect to r and introducing this expres- 
sion into Equation (12) with (6) and (8), we obtain a differential equation for 
A and r, which yields after integration 

A=l(l + (^) 2 ) ,r 2 :=a>0, (18) 
r 

where the boundary conditions are set such that A tends to when r — > oo. 
This relation is an extension of the expression (7) of A when thermodynamical 
fluctuations are taken into account. 

It is worth noting that this relation is strongly akin to a renormalization 
relation (e.g. see first Eq.(3.3) in [32], or Eq. (33) in [33], see also [34]). In 
the same way, Padmanabhan ([35], Eq.(19)) and Myung ([36], Eq. (14)) write 
a renormalized vacuum energy density by stressing a hierarchical structure of 
the form: 

^c = ^[l + (^) 2 + (^) 4 + ...], (19) 

where the first term is interpreted as the energy density necessitating to be 
renormalized, the second term as the vacuum fluctuations and the third one as 
the de Sitter thermal energy density (the dimensional constants and the factor 
3 were restored for convenience). Padmanabhan and Myung have introduced 
this form in the case L\ = Lp, here we use an intermediate scale L\ greater 
than Lp. A value of L\ accurate to observational results will be determined in 
Section 6. Setting Ai := j?, and multiplying the relation (19) by x, we obtain: 

A_-A 1 = -^[l + (^) 2 + ...]. (20) 



r| - r 2 
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Assuming Ai to be the bare CC Ao of the G.R. (eqs. (1) and (3)), this expression 
is of the same form as the expression (18) derived in our model. Let us discuss 
an intriguing point. The second term in the r.h.s. of (18) has by construction to 
be interpreted as a term of fluctuations, whereas it is the first one in the r.h.s. of 
(20) which is interpreted as fluctuations by Padmanabhan and Myung. Let us 
consider the solution to this apparent discrepancy. Relation (18) is derived from 
considerations concerning fluctuations by use of a Taylor expansion. It can be 
considered in the framework of rcnormalization as the truncature of a convergent 
series up to its first terms (see [33] after their Eq. (35)). In particular, Eq. (18) 
is valid for the two following limits: 

For r > ri, Eq (18) represents the first terms of a Taylor expansion and A 
tends to (7) which itself tends to its minimal constant limit A m = -|- when r 

r M 

tends to its maximum limit rM (IR cutoff). For these large values of r, we are 
concerned with the inflation undergone today by the universe. Therefore the 
fluctuations we considered in (14) and (18) are these recent fluctuations. In this 
limit, the model (7) was studied in [14]. 

For r<ri, Eq (18) can be rewritten as a Taylor expansion for little values 
of r 




and A tends towards another de Sitter state when r reaches its minimum limit 
r m (UV cutoff). For these values of r, what is concerned is the maximum 
value of A, during the inflation formerly undergone by the very early universe: 

Am = —ir- Padmanabhan and Myung consider these early fluctuations in (19). 

Considered as a renormalisation relation, Eq. (18) is also valid for any inter- 
mediate value of r and especially when r = r\. Assuming the given expression 
(18) for A, we derive now the evolution of the universe. 



5 Evolution equations for H and r 

In a spatially flat FLRW cosmology (k — 0, metrics signature H ), the field 

equations (2) for a two-component perfect fluid with the variable CC as first 
component are given by 

(A + X p)c 2 = 3H 2 (22) 
H=-l X c 2 P- (23) 

From now on, the index 2 for the second component of the fluid is understood. 
Equations (22) and (23) lead to 

H = -|(3i? 2 - Ac 2 ). (24) 

By defining 

Y := 1 , (25) 



8 



(24) and (10) become 

2H + 3-/H 2 - 3 7 c 2 r(l + aY) = (26) 

Y + 2HY - 2cY? = , (27) 

where the expression (18) of A has been introduced. Eqs (26) and (27) con- 
stitute a system of two coupled first order differential equations satisfied by 
{Y(t),H(t)}. We introduce the dimensionless variables 

Y H 

» := *'* := 17' (28) 

where Y\ and H\ are two, a priori independent, arbitrary constant scales for Y 
and H respectively. We also define the dimensionless time r 

r := (Hit)" , (n > 0). (29) 

Hence, (26) and (27) become 

2nT^ti + 3jh 2 -ZjXiyil + aiy) = (30) 

nr^y' + 2hy-2^/x' 1 yi =0 , (31) 

with 

Xi := (32) 
the dimensionless ratio of the scales Y\ and H\ and 

ai—aFi. (33) 
The apostrophe denotes the derivative with respect to r. Finally, we set 

A :_ * , ,34) 

where Ai is an arbitrary constant scale for A. With (18), (34) yields 

A=f(l + au/) (35) 

where 

is the dimensionless ratio of the scales Ai and Yi.We deduce 

U =Xlh (37) 
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which furnishes a scale for the dimensionless CC density parameter Da defined 

as 

Ac 2 

:= < 1 , (38) 

or, with (34)-(37) and (28), 

n A := Xl y -^yl. (39) 
Finally, setting r = Oa(0) = 1, with (39) gives 

M0) = [*iy(0)(l + aiy(0))]*. (40) 
Now it only remains to fix y(0), besides the parameters X\, a\, n , l\ and 7. 

6 Results 

We numerically solve the equations (30) and (31) and obtain the solutions y(r) 
and h(r), for the initial conditions (40) with 

t/(0) = 10 93 , 7 = 1, h = 1, X! = IO- 3 , ai = 10" 55 , and n = ^ . (41) 

By choosing r 2 = r 3 = rjj — 10 26 m (see Section 4, Eq. (17)), the condition 
(17) is satisfied for any r > r(0) ~ 10 _11 m (this last value will be derived in 
the following). 

The functions CI\(t) and A(r) are deduced from (39) and (35) respectively 
and the function Oa(t) is plotted in figures 1, 2 and 3 for different ranges of 
time. 

As the current values of the observational data (Qa(to) — 0.7 , H{tq) = 72 
kms~ 1 Mpc~ 1 ) of the KCDM model [37] [38] are model-dependent, we use the 
values derived in the framework of our K(t)CDM model from the studies of 
Snla observational constraints [14] 

ft A (r ) = 0.63 , H(t ) = 64 kms^Mpc' 1 , (42) 

where t represents the present instant. Let us emphasize that the use of the val- 
ues derived from the KCDM concordance model would only marginally modify 
the following results. 

The different scales H\,Yi, and Ai can be deduced from these two values 
(42). We deduce the instant r from the plot of the function Q(t) 

T = 0.9220725 . (43) 

Setting this value in y(r) and h(r) and using (35), we deduce 

3/(7-0) - 2.0576 x 10 18 , ft(r ) = 5.7149 x 10 7 , A(r ) = 2.0576 x 10 18 . (44) 
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Figure 1: Plot of the vacuum density parameter J7a as a continuous function 
of dimensionless parameter time t for r 6 [—0.930352, 0.930352]. The initial 
inflation (r = 0) and the final de Sitter stage appear. The function Q\(t) was 
parametrized by the interpolating functions y{r) and h(r), and the numerical 
solutions of Eq. (30) and (31) with the initial conditions (40) and (41). 



Then, Eq. (28) with the numerical values (42) and (44) yields 



H(t ) 
h(r ) ' 

From (32), (41) and (45), we have 



3.62928 x 10 



-26-1 



X H 2 

Y 1 = = 1.46555 x 1Q- J1 m- 2 . 



(45) 



(46) 



Besides, A(to) is directly fixed by the observational data (42) from (38) 

, . . 30a(t )H 2 (to) „ , 

A(r ) = — ^4 — — = 9-04649 x 10~ 53 m" 2 . 



(47) 



Thus with (34) and (44) 



Ai 



Afro) 
A(ro) 



4.39662 x 10~ 71 to 



(48) 



The numerical value of r\ linked to the scale Y\ and identified to the value L\ of 
Eq. (20) can now be estimated from (33), (41) and (46): r x = 8.26037 x 10 7 m. 

It is an intermediary scale between Lp ~ 8.02121 x 10~ 35 m (for l' p := \J^- — 

1.6 x 10~ 35 m, with L P = 2y/2nl'p ) and L Atoday = r(r ) ~ 1.82 1 04 x 10 26 m. 
We then obtain 



A(0) _ y(0)(l + a lV (Q)) 



A(t ) y(ro)(l + aiy(T )) 



4.86003 x 10 1 



(49) 
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Figure 2: Plot of Qa versus r at the end of the evolution. The universe enters 
again on a phase of acceleration for J7a > 1/3 and tends to de Sitter (plateau 
Oa = 1). The present time r = 0,9220725 i.e. to = 4, 88.10 17 s corresponds to 
^ao = 0.63 (see table 1). 

The dark energy density parameter fl\(r) decreases from the initial value f2A(0) = 
1 and then reaches the value 1/3, i.e. the end of the acceleration (or equiva- 
lcntly the exit of inflation), which occurs at the instant £2 = 1.23824 x 10~ 39 s 
(evaluated from Eq. (29)). 

Then, it passes through a minimum OA m (r m ) = 4.99970 x 10~ 20 at the 
instant t m = 2.9042 x 10" n s where A(r m ) = 8.795 x lO" 16 ™" 2 and H(r m ) = 
2.2955 x 10 10 s _1 . After, CI\(t) increases again ( A comes back [Krauss and 
Turner 1995 34b]) until a return of the inflation which occurs at the instant 
t 3 = 1.81 x 10 17 s = 5.73 Gyr . 

It reaches its present day value 0.63 at t (43) or t = 4.88084 x 10 17 s = 
15.4664Gyr. This value is in good agreement with the observational results 
derived from the "concordance model" ([38] Section 2.5. p. 3) and allows to 
avoid the age crisis (see e.g. [39]). With a linear extrapolation of H(t) during 
T3 — To, we roughly evaluate the redshift, namely z{t^) ~ 1.01826, where the 
acceleration begins again. This value is in good agreement up to the previous 
approximation with the exact value of the transition redshift z T = 0.965 found 
in [14]. 

In the same time, the "coincidence problem" is resolved. The date of the 
coincidence can be defined by different ways : for instance by the instant t m (see 
above), i.e. the instant from which f^A grows again, or by the instant tz (see 
above) when the universe inflates again, or even, more recently, by the instant 
t A = 3.61505 x 10 17 s = 11.4554Gyr when fl A = \ = fi m . In the following, we 
use this last definition. 

If we admit that the recombination starts at the instant i* = 1.17 x 10 13 s = 
372fcyr, and ends at = 1.54 x 10 13 s = 487kyr , the values of the main 
cosmological parameters are, for these instants : ^(t*) = 1.43xl0 -4 , £1\(t*) = 
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Figure 3: Plot of £1a versus r at the beginning of the evolution. There is an 
initial inflation (de Sitter = 1), then a graceful exit of inflation for < 1/3 
(see table 1). 



1.722 x 1CT 4 , H(t») = 5.68 x lO" 1 ^" 1 , #«) = 4.34 x Kr 14 s-\ A(r») = 
1.54 x 10~ 47 to~ 2 , A(t^) = 1.08 x lCT 47 m- 2 . A determination of the value 
of Qa^) ( or of Cl m {T*) ) from observational data, if feasible, would be an 
interesting test for the model. The following table presents the main results of 
the FLRW e.h. model. 
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T m = 0.6862 


tm = 2.90 lO" 11 


2.29 10 1U 


5 10-' 2U 


8.795 10~ ib 


r* = 0.8786 


t» = 1.17 10 ia 


5.68 10" i4 


1.43 10" 4 


1.542 10~ 47 


tI = 0.8796 


t'„ = 1.537 10 ia 


4.34 KT 14 


1.72 10" 4 


1.082 10~ 47 


t 3 = 0.9189 


t 3 = 2.29 10 iY 


3.52 10" i8 


0.333 


1.378 10~ b ' 2 


r 4 = 0.92081 


U = 2.87 10 iY 


3.15 10" 18 


0.5 


1.655 10~ 52 


t = 0.92207 


t = 4.88 10 iY 


2.07 10" 18 


0.63 


9.046 10~ 53 



Tablel. Values of H , fl\ and A at some key times of the evolution. 

The long period about from 6.39 x 10 _40 s to ~ 2.8 x 10 17 s when A does 
not dominate, i.e. when fl\ < 0.5, allows to maintain the standard model and 
its afferent results all along this period. 

Another interesting result to be noted is the shape of the curve CI\(t) ( 
fig.l) which presents a symmetry with respect to r = , suggesting a period 
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for Q\(t). Each curve ( £Ia, h, y, A, H, A) is (quasi-) symmetrical with respect 
to t = 0. The functions A(t) and £f(t) decrease (increase resp.) when r > 
(r < resp.). There is a "big bounce" on an inflationary space at t = 0, with 
a period that is stretched over more than twice to and another one at a very 
large time on a quasi- de Sitter space. We studied above the results given for 
r > only. The first inflation is "strong" (at t — 0) when A is maximum, the 
second one is "weak" (at t — > oo) when A is minimum (possibly^ 0). Today A 
is on the way to its minimum (maybe zero), and the universe is dominated by 
the small cosmological constant (37). 

We can also obtain the R(t) ( or "jerk" j(t) ) and S(t) "statefindcr" param- 
eters [40] and q(t) deceleration parameter for the model 

i + ^(i-n A ) (so) 



1 aH 2 2 
R=^ =3Xi yJ^^ { Vxj-H) + l (51) 



R 



(52) 



"3(g-|) ' 

with for their present day values 

qo = -0.445 , R = 0.610125, S = 0.137522 . (53) 

The positive sign of Rg indicates the present overinflation, (either increasing 
of inflation, or trend towards inflation). One can also evaluate R in the early 
period of deflation (i.e. desinflation) . For instance at r = 0.509275 when 
inflation was present, though decreasing, because Oa(t = 0.509275) = 0.4 > 
±, wc obtain R(t = 0.509275) = -1.414 < OWe can yet obtain R < for 
ft\ < 4 during some time after the end of the initial inflation. For example at 
r = 0.51151 when fl A = 0.2 < |, we have yet R(t = 0.51151) = -0.20005. The 
instant when R changes of sign is given by t 9 = 0.51204075 or tg = 3.07026 x 
10~ 39 s. Since then, R is positive. 

The comparison of the "kinematical" values (53) to observations (SNAP 
supernova data ) could also be a test for our model (see the remark of [40] , Sec- 
tion 4) and would allow the measurement of the production source of CDM [41] 
which is constrained by the statefindcr parameters and the "snap" parameter 

(4) 

* • off 1 ' 

Finally, this model also supports the low-^ CMBR power spectrum suppres- 
sion suspected in the WMAP data. Indeed, the event horizon plays naturally 
the role of a dynamical IR cut-off translatable to a wave number cutoff [42] [43] 
as for the holographic model [25] [30]. Applying the same procedure as [42] [43] 
for the Dirichlet boundary condition A c = 2r , we roughly evaluate the influence 
of this cutoff. In the integrated Sachs- Wolfe (ISW) effect, the IR cutoff appears 
as the lower limit of the comoving momentum which truncates the summation 
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giving the coefficient Ci of the power spectrum. For a power law spectrum with 
a spectral index n = 1 (Harrisson-Zeldovich), we have for the low I [42]: 

Ci^Nj2pf(kVo) (54) 

k>k c 

where ji are the spherical Bessel functions. In flat spacetime, the multipole 
number / is given by 

I = fcifab - V*) (55) 

where (?7o — f}*) is the comoving distance to the last scattering surface. The 
definition of the conformal time in flat space r\ — — - with the Dirichlet limits 
condition leads to 

i c = 7r((l + 0— -!)■ (56) 
Taking the current value z* = 1089 [5, 44] and using (18), we obtain 

r = ][^+f^)- (57) 
— can be evaluated with A*and A n from Table 1 and wc find : L = 5.15. For 

ro u c 

the value A'* of the end of the recombinaison, we get : l c — 6.76. These values 
are in the expected range. 



7 Conclusion 

From the FLRW c.h. entropy with thermal fluctuations (Eq. (14)), wc obtained 
the corresponding law for the cosmological constant A(r) as a function of the 
event horizon radius r (Eq. (18)). We used the striking similarity of this last 
expression (18) with a renormalization relation (Eq. (20)) to reinterpret (18) in 
terms of such a (truncated) renormalization relation (20). 

Using this A(r) law in an interacting two-component K(t)CDM cosmological 
model in a FLRW spacetime, we derived a system of coupled equations for H(t) 
and r(t) (Eqs. (30) and (31)). Through a numerical resolution (Section 6), we 
obtained the curves of evolution of J7a (figs. 1 to 3), the table 1 of the values of 
the main cosmological parameters {H, at different times in the evolution of 
the universe, and the evolution of the deceleration and statefinder parameters. 

The event horizon thermodynamical model presented here allows us to ac- 
count for most of essential features of the cosmology such as the value of 
the age of the universe, ~ 15.5Gyr (safe from the age problem [45]), the ra- 
tio ~5x 10 112 , as well as the exit of inflation after about t 2 ~ 23000 

Planck times, the existence of a minimum <~ 5 xl0~ 20 for Oa attained after 
about 3 xl0~ n s, the coming back of the inflation (after about 7.24Gyr from 
the beginning), followed by the recent coincidence, at 11.4Gyr, with the matter 
density, but also allows us to address the entropy problem, and the creation of 
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(dark) matter from the vacuum energy. Besides, the model does not perturb the 
results of the old standard model for the period going from 6.39 x 10~ 40 s~ 10000 
Planck times to 2.87 x 10 17 s = llAGyr all along which the DE is subdominant 
( f^A < 0.5 ). Finally, it can render an account of the low-? CMBR power 
spectrum suppression, for an angular scale corresponding to multipole numbers 
lower than l c ~ 7. 
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